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Abstract 

In this paper we show how to construct the coupled (multicomponent) Harry Dym (cHD) hierarchy 
from classical Stackel separable systems. Both nonlocal and purely differential parts of hierarchies 
are obtained. We also construct various classes of solutions of cHD hierarchy from solutions of 
corresponding Stackel systems. 
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0\ ■ 1 Introduction 

O . 

Various relations between finite- and infinite-dimensional nonlinear integrable systems have been inves- 
tigated since the middle of 70:s in a long sequence of papers starting from the paper pQ, through papers 
0-0 (see for example [3] for more detailed bibliography) and many others. In all these efforts, however, 
the main idea was to pass from infinite- to finite-dimensional integrable systems. This paper is a third 
paper in our series of papers showing that also an opposite way is possible: that of passing from ordinary 
differential equations integrable in the sense of Arnold-Liouville to infinite-dimensional integrable systems 
(soliton hierarchies). In paper [7J we demonstrated a way of generating commuting evolutionary flows 
from corresponding family of Stackel systems (that is classical finite dimensional Hamiltonian systems 
quadratic in momenta and separable in the sense of Hamilton- Jacobi theory) . We presented our idea in 
the setting of coupled (multicomponent) KdV hierarchies (for definition and properties of these hierar- 
chies, see for example [8]). In paper [9] we systematized and developed this idea by showing how solutions 
of these Stackel systems can be used for generating various classes of solutions of cKdV hierarchies. Al- 
though both papers have been written for the case of cKdV, similar constructions are possible for other 
hierarchies as well. In this paper we demonstrate a way of generating the coupled (i.e. multicomponent) 
Harry Dym (cHD) hierarchy (see [TO], [IT]) and various classes of its solutions from a class of Stackel 
systems of Benenti type. Our method leads both to the nonlocal cHD hierarchy as well as to purely 
differential cHD hierarchy, that is to a multicomponent generalization of HD hierarchy discussed in [12] 
(see also [IS])- The nonlocal part of cHD hierarchy has not been discussed in [TU] at all. We also clarify 
and simplify some of the results given in [7J, [5]. 

The paper is organized as follows. In Section 2 we briefly remind some basic fact about Stackel 
separable systems and discuss how they are related to corresponding Killing systems (dispersionless 
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nonlinear PDE's of evolutionary type defined by Killing tensors of Stackel systems). Sections 3 and 4 
are devoted to description of nonlocal multicomponent Harry Dym hierarchy and its various solutions, 
respectively. Sections 5 and 6 are devoted to local (purely differential) cHD hierarchy. 



2 Stackel systems and their dispersionless counterpart 

Stackel separable systems can be most conveniently obtained from an appropriate class of separation 
relations. Generally speaking, n equations of the form 

<#(Aj,jUj,ai,...,a n ) = 0, i=l,...,n, a.eR (1) 
(each involving only one pair A^,^ of canonical coordinates on a 2rt-dimensional Poisson manifold M) 
are called separation relations [14] provided that det ^ff 1 ^ 0. We can then locally resolve equations 
(fT]) with respect to <Zj obtaining 

di = Hi(\,(j,), i = l,...,n. (2) 

with some new functions (Hamiltonians) Hi{\, /1) that in turn generate n canonical Hamiltonian systems 
on M: 

X ti = —, »n = --^, . = l,...,n. (3) 

All the flows © mutually commute since the Hamiltonians Hi Poisson commute. Moreover, Hamilton- 
Jacobi equations for all the Hamiltonians Hi are separable in the (A, /i)-variables since they are alge- 
braically equivalent to the separation relations ([TJ. 

In this article we consider a special but important class of separation relations, namely 

n 

E a A" ' A <>- ' 7 A -- i = l,-,n (4) 
j=i 

with arbitrary fixed m, k € Z, e = ±1 (the constant j is not essential for the construction and is only 
introduced for a smoother identification our systems with the hierarchy in ( 10 )). The relations (U]) arc 
linear in the coefficients a* so that they can be (globally) solved by Cramer formulas, which yields 

<H = M T ^iG (m V + \v^ k) = H?' m ' k , i = l,...,n, m,keZ (5) 

where we denote A = (Ai, . . . , A„) T and fi = (/ii, . . . , fi n ) T . Functions Hi defined as the right hand sides 
of §5§ depend on m and k and can be interpreted as n quadratic in momenta /1 Hamiltonians on the phase 
space A4 — T*Q cotangent to a Riemannian manifold Q parametrized by (Ai, . . . , A n ) and equipped with 
the contravariant metric tensor G*-™- 1 (depending on m £ Z) given by: 

^-,...,^-J with A*= mXi-Xj). (6) 

It can be shown that G( m ) is of zero curvature for m = 0, . . . , n and that G^ n+1 ^ is of non-zero constant 
curvature, while all other choices of m lead to spaces of non-constant curvature. The Hamiltonians H^' m ' k 
are known in literature as Stackel Hamiltonians and the corresponding commuting Hamiltonian flows ([3]) 
are then called Stackel systems, or more precisely, Stackel systems of Benenti type. They are obviously 
separable in the sense of Hamilton- Jacobi theory since they by the very definition satisfy Stackel relations 
(0| . The objects Ki in ([SJ are Killing tensors for any metric and are given by 

( % % A 

K i = -A^y— j z = l,...,n, 

where qi — qi(X) are Viete polynomials (signed symmetric polynomials) in A: 

qi (\) = (-iy i = l,..-,n (7) 

l<si<S2<...<Si<n 
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that can also be considered as new coordinates on the Riemannian manifold Q (we will then refer to 
them as Viete coordinates). Notice that K{ do not depend on neither m nor k. Finally, the potentials 
can be constructed recursively [15] by 

vf +1) = vg\ - qi v} k \ k e Z, with V t {0) = S in , (8) 

(k) . . . (k) 

where we put V t = for i < or i > n. The first potentials are trivial: V t — (Ji.n-fc for k — 
0, 1, ...,n— 1. The first nontrivial potentials are = — For k > n the potentials become 
complicated polynomial functions of q. The recursion ([5]) can also be reversed 

V™ = V ^t 1] ~ —V£ k+1 \ k € Z, r = 1, . . . , n, (9) 

leading to potentials with k < 0. These potentials start with vf ^ = — an d are rather 

complicated rational functions of q. They will be referred to as negative potentials. It can also be shown 
[7] that 

g\f = yfto-m-i-j) (1Q) 

where = (G*-™* 1 ) 1 is the corresponding covariant metric tensor. 

Remark 1 The general n—time (simultaneous) solution for Hamilton equations f3|) associated with all 
the Hamiltonians |5P is given implicitly by 



-[= = r =d\ r , i = l,...,n. (11) 



U + a = ±- 

To see this it is enough to integrate the related Hamilton- Jacobi problem. Now, with n Hamiltonians 
j^n,m,k j n q we can assoc i a ^ ej by corresponding Legendre transforms, n Lagrangians £™' m ' TQ — >R 
given by 

L n l ^ k {\\ ti ) = -\lg^K~ 1 X u --X ik \ * = !,...,«. (12) 



Every Lagrangian £™' m ' fc give rise to n systems of Euler- Lagrange equations 



= 0, i = l,...,n (13) 



(each for every s between 1 and n) where 

d d 



Ej dXj dt s d{d\ 3 /dt s y J 

are components of the Euler-Lagrange operator with respect to the independent variable t a - 

Remark 2 By construction, the solutions ill]) are also general solutions for all the Euler-Lagrange equa- 
tions H13]) . It means that for a particular s the general solution of Euler-Lagrange equations Ej(L™' m ' k ) = 
is given by ill]) where t p for p ^ s plays a role of a constant parameter. 

Denote now the variable t\ as x (our method works similarly with any ti chosen as x). With every 
Killing tensor Ki for i = 2, . . . , n we can associate a dispersionless evolutionary PDE of the form 

A tj = KiX x = Zi [A] i = 2, . . . , n (14) 



(where A = (Ai, . . . A„) T ). We will call PDE's in (|14l) simply Killing systems. Here and in what follows we 
use the notation / [A] to denote integral-differential function of A i.e. a function of A, its x- derivatives and 
antiderivatives (integrals). In the case above Z\ [A] = Zi(X, X x ). The chosen variable t\ = x in (jT4")l plays 
thus the role of a space variable while the remaining variables U should then be considered as evolution 
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parameters (times). Equations ([Til) constitute a set of n— 1 integrable dispersionless equations that due to 
the form of Ki belong to the class of weakly nonlinear semi-Hamiltonian systems, i.e. hydrodynamic-type 
systems that are semi-Hamiltonian in the sense of Tsarev [TB] , [T7] and weakly nonlinear [TS] . Actually, 
the systems (|14[) are finite-component restrictions of the universal hydrodynamic hierarchy considered in 
[T5] . The variables A; are Riemann invariants of all the system (fT4"| as Ki are diagonal in A. The systems 
(fT4|) can also be considered as n — 1 dynamical systems on some infinite-dimensional function space V of 
vectors (Xi(x), . . . , X n (x)). with Zi being n — 1 vector fields on .M. It can be shown [18] that the vector 
fields Zi commute on V: 

[Zi,Zj]=0 i,j = 2,...,n. 

Note also that since K\ = I we can complete the system of equations (Tl"4")) by the equation A r = ifiA x = 
A K = ^with the translation-invariant general solution A^ = Xi(x + r). The vector field Z\ also commutes 
with all the vector fields Z%,. . . ,Z n [IS] . 

Proposition 3 Every mutual solution A(ii, . . . ,t n ) ill]) of all Hamiltonian systems (EP with Hamiltoni- 
ans of Benenti type (EJ) is (after replacing t\ with x) also a particular solution of all n — 1 corresponding 
Killing systems in |-?^| ). 

Proof. Let us assume that a vector function X(t%, . . . , t n ) solves ([TTj) . Then, by construction, it also solves 
the spatial part of ((3|) with appropriate functions (i(ti, . . . ,t n ) given by [a = dW(X, a)/dXi (W — W(X, a) 
is a common integral of all the Hamilton- Jacobi equations for Hamiltonians ff™> m ' fe ) jt means that 
X(ti, . . . , t n ) solves 

X u = -^H? m ' k = ZKiGWn, i = 1, n. (15) 

Since Ki = I we get from the first equation in (fT5|) fi(ti, . . . , t n ) = hg( m '(\(ti, . . . ,t n ))X tl (ti, . . . , t n ). 
Substituting it to the remaining equations in (|15l) yields then 

X u (h 7 . . . ,t n ) = Ki{X{t u . ..,t n ))X tl {ti,...,t n ), i = 2,...,n 

which concludes the proof as t\ — x. Thus, all the solutions (jTJJ) also solve all n — 1 Killing systems in 

(UJ). ■ 

Moreover, we have 

Theorem 4 The general (n-time) solution of all the Killing systems in \1J$ is given by 

ti + Ci = y^l — j—dX r , i = l,...,n (16) 

f^J <^r(A r ) 

( where (p r are arbitrary functions of one variable ) 

The proof of this statement can be found in [T8j . Obviously, (Tf6|) contains all the solutions (fTT|) . 

Suppose now that a particular solution (|16j) of our Killing systems ()14() is of the more specific form 
pTjl . Since this class of solutions - by construction - satisfies all the Euler-Lagrange equations (TT5)) . we 
can treat equations (JT3J) as additional bonds that these solutions satisfy. We can therefore use these 
bonds to express some variables A^ by other A's. Thus, within the class (TTT[) of solutions flTBl of Killing 
systems (|14[) we can perform a variable elimination (reparametrization) that turns (|14l) into entirely new 
sets of evolutionary PDE's. As we have demonstrated in [7] and in [5], in carefully chosen cases and 
in a particular coordinate system (Viete coordinates ([7)1) this reparametrization turns systems (| 14[) into 
systems with dispersion (soliton hierarchies) with the solution (1111) being also a solution of these new 
systems with dispersion. In this paper we will produce by this method (the local and the nonlocal part 
of) the coupled (multicomponent) Harry Dym hierarchy. 
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3 Nonlocal coupled Harry Dym hierarchy 



Assume now that e = 1 in (j4|) and therefore also in ([TTj) . (|12|) etc.). In order to perform the elimination 
procedure just mentioned, let us pass to Viete coordinates as given in ([7]). The Killing systems (fbl)) are 
tensorial so in Viete coordinates they have the form 

Qu = Ki(q)q x , i = 2,...,n 

or, explicitly 

d j ~ 1 

■jrQj = (<b+i-i)x + ( qk (lj+i-k-i) x - Qj+i-k-i (qk) x ) = (Z? [q\) J , j = l,...,n (17) 
1 fc=i 

(where we put q a = for a > n), where i = 2, . . . , n and where (Z™ denotes the j-th component of 
the vector field Zi [q] . The superscript n at Z^ indicates the number of components in the vector field Zi 
and we will sometimes use it since we will need to switch between various n. From (|17|) one can see that 
{Zf [q]Y = (ZJ [q]Y for all i,j = l,...,n. Obviously, G( m \ and Ki are tensors and can thus also 
easily be transformed to Viete coordinates. 

Consider now Euler-Lagrange equations (TT51) with s = 1 (so that t s = t\ = x) associated with 
Lagrangians ^• m ' k denoted further on for simplicity as L n ' m - k . Denote also E\ as J?j, i = 1, . . . ,n and 
consider the equations 

E i (L n ' m ' k )=0,i = l,...,n, n,eN, fceZ, (18) 
written in q- variables, so that now 

Eji — „ , „ , ? — i, . . . , n, 

while (since K\ = I) 

L n, m .k = L n>m,k^ qx) = l q T g (m) qx _ ^W. (19) 

As it has been shown in [7] the following symmetry relations are satisfied for a = 1, . . . , n — 1 

( L n,m,*) = E ._ a ( L n,m+a,k-a-) ) 4 = a + 1, n , (20) 

that can also be written as 

Ei (L n > m ' k ) = E l+a (L n < m - a < k+a ) , i = l,...,n-a. (21) 

Due to (|20p and (|2"Tj) the equations (TT51) can be embedded in the following double-infinite multi-Lagrangian 
" ladder" of Euler-Lagrange equations of the form 

Ex = E2 ( L n,m+i-2,fe-j+2^ = . . . = j^n.m+i-n.fc-i+n) = Q) j = . . . , -1, 0, 1, . . . 

(22) 

with fixed m, fc £ Z (the equations (fT5|l fit in (|2"21 at j = 1, 2, . . . , n). For a given dimension n the ladder 
(|22p is determined by the sum m + k in the sense that various choices of m and fc with the same m + k 
yield the same ladder. 

We are now ready to present our elimination procedure leading to multicomponent integral (nonlocal) 
Harry Dym hierarchy. Assume that we want to produce first s — 1 flows of the TV-component (N £ N) 
hierarchy. Let us take n = s + N — 1, m = — N and k = in (fT2")l . that is, let us consider the purely 
kinetic Lagrangian £ J n -" N - Q with n = s + N — 1 and the corresponding Euler-Lagrange equations (fTE)) . 
Due to this special choice of all parameters the last n — N equations in (I18[) attain the form 

E N+ i (L n - N >°) = -jq, hxx + <p n -N[qi, —,Qn-l] = 0, 

: " (23) 

B„_i (I n " w '°) = + V2[?l, <7N+l] = 

£ n (L"--^ ) = -igjv+i,xx + Vi[gi,-,9Jv] = 0. 
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and are a part of the ladder (f2"2l with m + k = —N. Now, by direct calculation of Ei (L n > '°j with the 

(i) 

use of some identities satisfied by the potentials V\ it can be proved that 

P (rn,-N,0\ _ ei ( T n+l,-N,0\ , 1 

i'iV (-^ J — {-Li ) + -q n +l,xx, 

E . = (in+i.-jv.o) ; i = jV + 1, . . . ,n. 

These identities lead to 

Proposition 5 The functions tfi in i23\) do not depend on n in the sense that increasing n to n + 1 (and 
keeping N constant) turn f23\) into n — N + 1 equations 

En+i (L n+1 '~ N '°) = — j<ln+i,xx + En (L n ~ N '°) = — \q n +x,xx + fn-N+ilqi, In] =0, 
E N +2 (L"+ 1 -- Ar -°) = E N+1 (L n >- N >°) = -lq n ,xx + Vn-N[qi,-,qn-l] = 0, 

: : : (24) 

E n (L n+1 '- N >°) = £7 n _i (L n '~ N '°) = -hqN+2,** + <P2[qi,...,qN+i] -0, 

E n+1 (L"+i.-w.«>) = £ n (i™.-^°) = -|^+i ira + ^i[gi,...,9jv] -0. 

It means that increasing n to n + 1 (and keeping TV constant) in (f2"5)l dcs not alter these equations 
except that a new equation of the form 

£at+i (£™+i,-.iV,o^ ^ __ qn+1 ^ x + ( ^ n _ JV+1 [g 1) ...,g n ] = 

is added at the top of (|2^|) . As we will see soon, this will result in the fact that our construction indeed 
yields an infinite hierarchy of commuting flows. 

Due to their structure, equations (|23p can be formally solved with respect to the variables qN+i, • • • > 9n; 
which yields qN+i, • • ■ > <7ti as some nonlocal (integral-differential) functions of qi, . . . , q^' 

qN+l = fl [qi,...,q N ] 

: (25) 
q n = fn-N+l [qi, ■ ■ ■ , qN] , 

where, due to Proposition [SJ the functions do not depend on n, so increasing n by 1 (and keeping 
N constant) will only result in one new equation at the bottom place in (|25[) . Let us now replace the 
variables qN+i, ■ ■ ■ » qn m the first N components of the first s — 1 Killing systems (|17|) by the corresponding 
functions /, (right-hand sides of (|2"5jl ). This yields equations of the form 

N 

q tr =Z r [q] r = 2,...s (26) 

where q denotes the first N entries in q i.e. q = (qi, . . . ,qN) T ■ They are in general highly nonlinear 
autonomous systems of N evolution equations for q\, . . . , qN- 

Theorem 6 The vector fields Z r [q] in S26\) do not depend on s in the sense that if we increase s by one 
in our procedure then \2b}) are unaltered and a new equation q ts+1 = Z s +i [?] appears. 

Proof. This theorem is a consequence of Proposition [5] If we increase s to s + 1 and keep TV constant we 
have to take n + 1 instead of n in our procedure as n = s + N — 1. Due to (117|) we have (Z? +1 Y = (Z?) j 
for r = 2, . . . , s and for j = 1, . . . , N i.e. the first N components of the first s — 1 of Killing systems (fTT]) 
do not change when we increase n to n + 1. Moreover, as we explained above, the n — N functions f t in 
(f2"S"]) do not change either. So, the elimination procedure for the first s — 1 vector fields Zi is not altered 

N _ ,1 

leading to exactly the same vector fields Z r [q] with r = 2, . . . , s while the vector field Z££i yields the 

vector field Z^ +1 [g] i.e. a new equation at the end of the sequence ([26)) . ■ 

Repeating this argument we can increase s indefinitely. Thus, our procedure leads to an infinite 
hierarchy of evolutionary vector fields (flows) 

q tr = Zr[q] r = 2, 3, . . . (27) 
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in the sense that if we wish to produce any first s — 1 flows (|26[) of the hierarchy we can perform our 
procedure with n = s + N — 1. This way we can obtain arbitrary long sequences of the same infinite set 
of vector fields with dispersion that pairwise commute (soliton hierarchy): 

Theorem 7 The vector fields Z^ \q] commute i. e. 



for any i,j = 2,3, . . 



This theorem is due to the fact that the original vector fields Zf commute and that the Euler-Lagrangc 
equations Ei{L^ n ' m,k ^) — are invariant with respect to all the fields Zf [7]. Moreover, the vector fields 

Z i still commute with Z 1 = (qi tX , . . . , Qn,x) ■ As we demonstrate below, the hierarchy (|27|) is the 
nonlocal part of the multicomponent Harry Dym soliton hierarchy as discussed in |12) . 

Example 8 Consider first N = 1 (one- component hierarchy as discussed in [12]). Suppose that we want 
to obtain the first s — 1 = 2 flows of the hierarchy. We have then to take n = s + N — 1=3 and consider 
the elimination equations \23\) for these parameters. The pure kinetic Lagrangian L 3 '~ 1,0 has the form 

L 3 = ^qj x \q x q 2 - ^93 - ^qf \ + ~qi, x Q2,x (<3i - 92) - j^i.^i - ^l^.x^i + ^?2,x9V 
so that h!3\) become 

e 2 (l 3 ' -1 ' ) = ~q 3 , xx + ^qi, xx q2 + ~92,xx9i - ~9i,xx9i - T^qiqlx + \fl\,x(h,x = o, (28) 

£ 3 (L 3 '- 1 ' ) = ~q 2 , xx + ~9i,xx9l + ~ql >x = 0. 

Due to their specific structure, we can solve &28\) with respect to q 2 and q 3 . We will thus use i28\) to 
eliminate variables in the corresponding n = 3-component Killing systems |J7[ ) that have in this case the 
form: 

92,x 

Q3,x + qiQ2, x - 9291.x I = z$, 
9193,2 - 9391,2 







(:) 


-1 










(:! 









q3,x 

x - I 

9293,2 - q 3 q2 



qiq3,x - 539i,x 1 = z\ . (29) 



By the second equation in (28\) we obtain 



1 2 U_ 

92,2:2; — 2^1,2; + qi,xxQi- 



1 L Vl 

92,2; = ^qiqi, x + 2° 9i 



d- 1 = / ...dx + <p(t 2 ,t 3 ) 



Integrating it once we obtain 



where 



is the integration operator with the integration parameter if that has to be chosen from case to case and 
has therefore to be treated as a part of the solution of every integration problem. It is always possible to 
find such a function. Integrating q 2 x we obtain 

1 2 1,-2 
92 = ^9i + 9i9i,xx- 

Further, the first equation in (28\) yields 

93,2-2; = 91,2-2:92 + 92,2:2:91 - 9i, 22-9i _ 9i9i x + 91,292,2 (30) 
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Inserting to it q 2 and q 2iX , calculated above, and integrating once we obtain 

<l3,x = —^d~ l qfqi tXX + ^qid^ 1 q 1 q ltXX + ^q 2 qi, x + ^qx,xd~ 2 qiqi,xx (31) 

By inserting the obtained formulas for q 2tX and q%^ x into the first N = 1 components of Z 2 and Z% we 
obtain the first two flows of our nonlocal soliton hierarchy: 

<li,t 2 = ^qm, x + ^d^qxqi^x = Z 2 , (32) 

qi,t 3 = -^d~ 1 qlq 1 , xx + ^qid~ 1 q 1 q 1 . xx + -^q\qx,x + ]^qi, x d~ 2 qiqi, xx = Z 3 . 

Observe that in this particular case we did not have to calculate q 3 since it does not enter into the first 
component of neither Z 2 nor Z 3 . We needed howcvev Q2 27i order to calculate Qs,x • 

The flows \EM 

commute due to Theorem^ 

Example 9 Let us now take N — 2 and s — 1 = 1 so that n — 3 again. We will thus eliminate n — N = 1 
variables ( namely q 3 ) from the first N — 2 components of the field Z 2 above. The elimination equations 
(03) reduce now to E 3 (L 3 ^ 2 > a ) =0. But, according to $2$, E 3 (i 3 ^ 2 >°) = E 2 (L 3 ' -1 ' -1 ) = E 2 (L 3 ^ ) , 

the last equality due to the fact that L 3 ^ 1 ^ 1 = L 3 '~ — = L 3 '~ 1,0 + — Thus, the elimination 

equation E 3 (L 3,_2 '°) = coincides with the first equation in S28\) and yields exactly S30\) . Plugging its 
integrated form h31\) into the first two components of Z\ yields the first flow of the 2-component nonlocal 
cHD hierarchy: 



d 



q\ \ I qi, x 



The map 



(33) 



Ul = E N _ i+1 (L N >°>°) , i = l,-..,JV (34) 

transforms the hierarchy (|2"T)l into the nonlocal part of the coupled Harry Dym hierarchy (see [TU] for 
its local part) that is the generalization of the one- field nonlocal HD hierarchy presented in [12"] . For 
example, for N = 2 this map reads 

1 

u i — ——qi,xx 

1 1 2 1 

u 2 = --^qi, X x + -^qi, x + -^qiqi,xx 

and applied to the field Z 2 above yields 



( Ul \ — of u i,xd 2 ui + 2uid 1 m - \u 2: .. 



= —2 

dt 2 \ u 2 J y u 2 ^ x d~ 2 ui + 2u 2 d~ 1 u 

4 Solutions of the multicomponent nonlocal HD hierarchy 

We will now construct a variety of solutions of the hierarchy (|27|) . 

Theorem 10 For any fj e {0, . . . , n — 1} the functions Xi = A,(ti, . . . , t n ) given implicitly by 



1 _ r \ n ~ l 
= ±-J2 —= r d\ r , i = l,...,n. (35) 



' 1 i/A, i^.[ =1 aj\ r J -±\ 

are such that the corresponding functions qi — qi(x — t±, t 2l . . . , t n ), i — 1,...,N, given by ^ are 
solutions of the first n — fj (n — 1 for fj = 0, 1) equations of the N -component integral cHD hierarchy 
\21^ . The variables t 2 , ■ ■ ■ , t n —p+i (t2,---,t n for fj = 0,1) play then the role of evolution parameters 
(dynamical times) while the remaining U 's are free parameters. 
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For the proof of this theorem, see Appendix. We will now consider some particular, interesting classes 
of solutions (|3"5)) . Assume that /3 = in (|3"5"j) and that aj — j5j. n + <5j,n-7 for some 7 € {0, . . . , n — 1}. 
Then ([55)) attain the form 

t l + c t = ±-J2j^==dX r , i=l,...,n, 

that integrated yields 

1 n 

f. + r . —4- >p \ n-i+N/2-7/2+1 ■ 1 /o fi \ 

The above system can be algebraically solved with respect to Ai only for two choices of 7, namely 7 = N 
and 7 = A + 1, but it turns out that the case 7 = A leads to trivial solutions (polynomial solutions not 
depending on x). Thus, we must assume 7 = A + 1. In this case the above equations attain the form 

n 

^ = ± 2(n-i + l/2) £ A "~ W/2 > i = 1 >-' * < 37 > 

Note that ([37]) do not depend on A. It means that for any A between 1 and n — 2 (as7 = A + l < 
n — l)the functions q±(x,t2, . . ■ , t n ), ■ • ■ , qN(x, *2, . . . , t n ) obtained from ([37| through ([7]) solve the first 
n — 1 equations in (|27[) . The following two examples illustrate this. 

Example 11 Assume that n = 3. Then \3ty attain the form (with x = *i, Ci = 0, we also choose only 

+ m gg); 

X = ^T,i=l Z i = \ (Pl ~ 5 (PlP2 ~ Pz){p\ - P2)) 

h = \Y*=x4 = l(pl-3 PlP2 + 3p 3 ) (38) 
*3 = ELi z » = Pi 

where Zi — A^ 2 , i = 1,2,3 and where p\ = Y^i=i z i> P2 = ^1^2 + 2i#3 + Z2Z3 and p 3 = z\z 2 z 3 are 
elementary symmetric polynomials in Zi . The right hand sides of LS8\) follow from Newton formulas: 



E?=i*r= E ( _ ir+ ^e+- m K + °* + • • • + «n 1)1 ?1 a a p „ n /or ro < n , 

ai! . . . a n ! 

Qi +2a2 +. . — 771 

(39) 

expressing sums of powers of variables as functions of their symmetric polynomials ( these formulas can 
easily be extended to the case m > n by taking n' — m and putting all p n +i, • • • , p m equal to zero). The 
system can be solved explicitly yielding the solution flff7| ) in p-variables: 



Pi = *3 

-T5z-2*j+15*j*2 
1,1 5(3* 2 -*|) M,n 

15t 2 *l + 45t 2 - *3 - 45x* 3 

P3 = 



15 (3*2 - *!) 

On the other hand, according with and with \39\) 

Qi = -(Ai + A 2 + A 3 ) = - (z 2 + z 2 2 + 4) = -{2p 2 - p\). 
Plugging [JW into the above identity we obtain 



4 + 15*|*2 - 30a; 
5(3*2 -*!, 



?i(x,*2,t 3 ) = qi{p l (x,t 2 ,t 3 )) = 3 -j— j 73T . (41) 



According to Theorem I * 01 the function Qi(x, *2, *3) given by |^*[ ) yie/d a two-time solution to the first 
n — 1 = 2 flows of the nonlocal 1-field (i.e. with N = 1) HD hierarchy i.e. to both systems 

(after an appropriate choice of integration constants). 
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Example 12 Let us now take n = 4. In this case the equations JffTp read (again wit all Ci — and 
+ only and due to i39\) ) 

x = ^Y,t=i z I = -j{p[- 7(PiP2 ~ Ps) ((pf - P2) 2 + pi/33) - 7p 4 (p? - 2pip 2 + ps)) 

*a = gEt=i*< = J (Pi - 5 (PiP2 ~ Ps)(p? - P2) - 5 PlP4 ) (42) 

*3 = 3X^=1^ = 3 (Pi - 3piP2 + 3p 3 ) 

*4 = Ei=l^ = Pi 

where as before Zj, — X^ 2 and pi are again symmetric polynomials of the variables z%, . . . , Z4. This system 
can again be algebraically solved for p\,...,PA although the solutions are too complicated to present them 
here. We have now, according with 

gi = -(Ai + A 2 + A 3 + A 4 ) = - {z\ + zj + z% + zf) = -(2/92 - pj) 
q 2 = AiA 2 H h A3A4 = z\z\ H h z\z\ = p\ — 2pip 3 + 2p 4 

Substituting the variables pi obtained by solving into these expressions we obtain expressions for 

qi(x,t2,t 3 ,ti) and qi(x,t 2 ,h,t 4: ): 

I . . . x Pl(x,t 2 ,t 3 ,t 4 ) P2(x,t 2l t 3l t 4 ) 

qi{x,t 2 ,t 3 ,t 4 ) = , q 2 (x,t 2 ,t 3 ,t 4 ) = n2( , - - x (43) 

<2(*2,*3,*4j Q Z {t 2 ,t 3 ,ti) 

where Pi and Q are rather complicated, but perfectly manageable for any computer algebra program, 
polynomials. More specifically 

P 1 {x,t 2 ,t 3 ,t 4 ) = ~ (105^2 - t! - 21*1*3 + 630*2*3 - 630xt 4 + 315*§*|) 

and 

Q(ta, h, U) = 45t 2 U +t%- 15*3*1 - 454 

while P 2 is a quadratic in x polynom that is too complicated to present it here. Now, according to Theorem 
[77J1 and the theory above, the function q\{x, t 2 , t 3 , t 4 ) in \43\ ) solves the first n — 1 = 3 1- field flows of the 
hierarchy {27\) , so in particular both the flows 132\), while the vector function 

qi{x,t 2l t 3l t 4 ) 
qi{x,t 2 ,t 3l t 4 ) 

solves the first n — 1 = 3 flows of the N — 2-field cHD hierarchy starting with V33\) . 

Let us also remark that formulas (|3T>1) often lead to implicit solutions of ([77]) . We illustrate it in the 
following example. Choose N = 1, n = 2 and 7 = in (|36|) . This yields (again for Cj = 0) 

x = ~ (4 + zl) , t 2 = i (zf + zl) (44) 

(with Zi = X 1 / 2 ) that can not be algebraically solved. However, (|44"|) can be embedded in the algebraically 
solvable system (|3"51) in the sense that (f3"5f reduces to (|4"H) if we put z 3 — or equivalent ly p 3 = 0, since 
p 3 = z~iz 2 z 3 . By virtue of Theorem 1101 it means that the function 

j/(M 2 ) 5 + 152/(:M2) 3 t2 - 30x 

q 1 {x,t 2 ,y{x,t2)) = rTTT 7 , x 3 x 

5 (3i 2 - y{x,t 2 y) 

with the variable y(x,t 2 ) defined implicitly by the equation 

15t 2 y 3 + 45*| - y 6 - 45xy = 

(i.e. by the last equation in (|4"0")) with y instead of t 3 ), also satisfies the first flow of the nonlocal HD- 
hierarchy i.e. the first flow in 
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5 Differential (local) cHD hierarchy and its solutions 



We will now obtain the purely differential part of cHD hierarchy as well as a class of its implicit solutions. 
We choose now e = — 1 in (0} in order to obtain real solutions in the local case (note that it does not 
influence the potentials ). Analogously to the case of nonlocal hierarchy, we will perform some 
variable elimination on the sequence of Killing systems (1171) . Suppose thus that we want to produce the 
first s flows of the ./V-componcnt local (i.e. purely differential) Harry-Dym hierarchy. Put n = s + N and 
consider the first n — N Euler-Lagrange equations for the Lagrangian £™<™ _Ar ' _n . Using the fact that 
Vy = Vy (q n -j+i, . . . ,q n ) it can be shown that they attain the form 

Ei (L n > n - N >- n ) ee -L + 7 W [ qi ,...,q N }= 0, (45) 
Ei (L n < n - N >- n ) ee + 7 W . . . , qN _ i+l] + _l_ 7 g ) [qn _ %+2: ...,q n }=0, i = 2, n-N. 

qn 

where as usual q a = for a < 1. Note that P5|) and ([25)1 belong to the same ladder (p?Hj) of Euler-Lagrange 
equations since in both cases m + k = —N. 

Proposition 13 The functions 7 ji )7 ,j2 , 7l c?o ^oZ depend on n in the sense that increasing n to 
n + 1 u/iZZ not aZter J^5p except that a new equation originates at the bottom of the list |^5[). 



The proof of this proposition resembles the proof of the analogous statement for nonlocal case i.e. 
Proposition [5j Note now that the structure of (l4"5"j) makes it possible to eliminate (express) the variables 
ojv+i, . . . , q n as (purely differential now) functions of o/i, . . . , qjy (although now, opposite to the nonlocal 
case, we first calculate q n , then q n -x and so on up to q^+i)'- 

Qn = f[ N) tol,---,<?iv], 

: (46) 
qN+i = fn [qi,...,q N ]- 

Now, let us replace the variables qN+i, in the first N components of the last s systems in (|17p. 

That leads to s highly nonlinear (purely differential) evolutionary equations of the form 

q tr =Z^[q] r = n-s + l =N+l,...n (47) 

where as before q = (?i,-.-,?jv) but with new, purely differential, vector fields Z r . These fields 
constitute in fact the first s fields of the local cHD hierarchy. Contrary to the nonlocal case, however, 
the first field of the hierarchy appears as the last equation in (|4T|) i.e. q t = Z n [q], the second field is 

q t — Z^-i [q] and so on so that the fields of the hierarchy originate in (|47|) in the reverse order. We 
will therefore introduce a new notation and denote 

t p = t n -p + i, and Xp = Z„_ p+1 , p = 1, . . . ,n - 1 (48) 
so that q tn — Z^ [q] reads q Ti = Xy [q] and so on. The sequence ((47]) becomes therefore 

TV 

q Tr =X r [q], r = l,...,s. (49) 

A theorem analogous to Theorem [6] explains that this procedure leads to a hierarchy. 

Theorem 14 The vector fields in do not depend on s in the sense that if we increase s to s + 1 
then the above elimination procedure produces the same sequence |^ff[ ) of evolutionary systems plus a new 
system q Ta+1 = X s+1 \q] at the end of the sequence \4ty (i- e - °t the beginning of the sequence j47\)). 
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Proof. Consider the s systems (jTf)) and increase s to s + 1 keeping N constant. We have then to 
take n + 1 instead of n in our elimination procedure. Since, according to Proposition 1131 the functions 
7i l ' TiaP 1 7i do not depend on n the functions do not depend on n either. It means that increasing 
n to n + 1 (and keeping AT constant) turns the equations (|46p into 

9n+l = /i [9i,---,9at] 

= fn N) [qi, ■ ■ ■ , Ojv] 
QN+1 = fn+\ [qi,-- -,Qn] 

and at the same time the the structure of the last s equations in (fTTf changes so that q n is replaced by 
qn+ii 9n-i is replaced by q n and so on until qN+2- It means that the last s equations in the (extended to 
n + 1) sequence (|47[) will after elimination remain the same while a new equations originates - this time 
before (with lowest r) the other s ones. ■ 

Thus, by taking appropriate s we can produce on demand an arbitrary (finite) number of evolutionary 
vector fields 

q T =X r [q], r = l,2... (50) 



and due to same argument as in the nonlocal case, these vector fields all mutually commute: 

= for all i,j = 1,2... 
The described procedure leads in fact to multicomponent local Harry Dym hierarchy. 



Y N -Y N 

i ) j 



Example 15 Let us first produce the first s — 2 flows of the standard Harry Dym hierarchy i.e. with 
N = 1. We have n = s + N = 3. Consider the Lagrangian 

T n,n-N -n _ r 3,2.-3 _ * 2 %x^x 93^x92 qi _ q 2 

and the corresponding Euler- Lagrange equations J^5[ ). TTiey attain the form 

Bk - ^ - ^9i,xx - 

1 j 2<z 3 3 4<z 3 2 ^ 2 93 

and can i/ms easz/y ae solved with respect to q2 and q% yielding ffllfy in the explicit form 



XX ) 



.-1/2 



93 = 93 [qi] = (2gi. 

92 = q%[q\] = 5 (59l,xxa! - 4 9l,xx9l,xxxx) (2gi,a;a) _7/2 

Substituting these expressions to the first ( since N = 1 ) component of the last s = 2 Killing systems of 
the sequence |17| ) we obtain the following two commuting flows: 

9i,*2 = (92[9i])x > 9i,t 3 = fefeiDz 

or 

9i,T! = (93[9i]) a = , 9i,r 2 = (9a[9i])x = X 2 

(with the differential functions 52(91] and q^ [q±] given as above) i.e. the first two members of the well 
known local Harry Dym hierarchy. 
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Example 16 Let us now produce the first s = 2 flows of the N — 2- component Harry Dym hierarchy, we 
need therefore n = s+N = 4. The Euler- Lagrange equations |^5[ ) for the Lagrangian L n i n ~ N <~ n = L 4 - 2 ~ 4 
attain the form 



^4,2.-4x1,1,^ 1 2 J 



2oj 2 y " 



ttaf is soluble with respect to qs and q^ yielding 

qi = g4[(7i,(?2] = -w _1//2 = - (2q 2/X x - ql iX - 2qiqi,xx) 1/2 

q3 = qi[qi,q2\ = ~qi : x X w~ 3/2 

Substituting these functions to the first N — 2 components of the last s = 2 Killing systems of the sequence 
|i7[ ) ( with n = A) yields the desired flows 

and 

d (qi\_x 2 = { (qi,xxW- 3 / 2 ) x \ 

dr 2 \q 2 J~ 2 ~\ :i V. - 9i + (w^/ 2 ) x ) ^ 

Our parametrization of Harry Dym hierarchy differs from the parametrization given in |10j . Generally 
speaking, the hierarchy (|50[) is transformed into the multicomponent Harry Dym hierarchy presented in 
[TU] through a complex version of the map (|3"4")l 

u r = -iSjv-r+i (L N '°-°) , r = 1, • • • ,N, i 2 = -1. (53) 
For example, in the u-variables the system (|5ip attains the form 

d ( U ^\-Y 2 ^=( '("^Lx 

*T U J ' 1 M ' ^ ux (V 1/2 ) r + K l/s V 

that is exactly the flow (24a) in [TU]. 

We will now formulate a theorem corresponding to Theorem I1Q[ i.e. we will generate a wide class of 
solutions of the hierarchy ([50]) . 



Theorem 17 for any ft £ {0, . . . , n — 1} the functions Xi = A»(ti, . . . , t n ) given implicitly by 

^ + c i = ±-V / r dA r> i = l,...,n. (54) 

are such that the corresponding functions qi — q%(x — t\, t 2 , ■ ■ ■ , t n ), i = l,...,N, given by |7p are 
solutions of the first n — ft (n — 1 for ft = 0, 1) equations of the N -component integral cHD hierarchy 150\) . 
The variables t/3+1 = T n —p, . . . ,t n = T\ (t 2 ,...,t n for ft = 0,1) are evolution parameters (dynamical 
times) while the remaining t% 's are free parameters. 

We will not prove this theorem here as its proof resembles the proof of Theorem [TDJ Comparing 
Theorems [TOj and [17] we can see that the solutions ([35]) and ([54]) are for ft = 1, . . . , n — 1 related through 
the transformation ft — > n — ft, e — > — e. i.e. every solution ([35]) for ft — l,...,n — 1 coincides, 
after changing e — > —e, with the solution ([M]) with ft' = n — ft. It also means that the nonlocal flow 

qt n _p +1 = Zn-p+i [9] an( i the local flow q Tfj = Xp [5] share the same family of solutions, namely (133]) (or 
([51]) with ft' = n - ft and with e' = -e). 

It turns out that (|54|) cannot be explicitly solved. However, by taking all a, = in ([54]) (which yields 
the so called zero-energy solutions) we can obtain interesting implicit solutions to our hierarchy ([50]) . 
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Example 18 Consider the solutions \5J$ with N — 2, n = 3 and with all a, = 0. They have the form 

3 

U+ Cl = ±^2 X 4 - l d\ r , i = 1,2,3. (55) 

r=l ^ 

fi/ie same /or /3 since j3 -terms cancel after inserting ai — 0) and according to Theorem \ 1 7| they solve 
the first n — 1 = 2 flows of the N — 2- component cHD hierarchy i50\) i.e. both the flows k51\) and \52\) . 
Equations \55)) after integrating yield (remember that t\ = x; we also put all Cj = for simplicity of the 
formulas ) 

z=j£L^ h = ^ELi X l *s = \ELi X l (56) 

which can not be algebraically solved. However, similarly as in the nonlocal case, we can embed 156]) in 
the system 

h = ^EtiA? = I (pf - 3pip 2 + 3 P3 ) (57) 
h = Z)i=i^< = Pi 

(where pi are symmetric polynomials in Xi so that = (—1)* Pi) in the sense that putting A4 = (so that 
P4 = since p^ — A1A2A3A4; the righ hand sides of |57| ) are again due to iS9\) ) in J<57| ) we obtain 156]) . 
The equations J57| ) can be explicitly solved yielding. 

qi = -pi = -U 

1 2 

92 = P2 = ~*3 + 

1 3 

93 = ~P3 = -*2 - g*4 + *3*4 

1 4 1 2 1 

94 = P4 = -X + —t 4 - -t 4 t 3 + -t 



t 2 t4 



Thus, the functions qi(x,t2,t 3 ) given implicitly by 

qi = -pi = -t 4 (x,t 1 ,t 2 ) 

92 = P2 = -t 3 (x,t 1 ,t 2 ) + -t4,(x,h,t 2 ) 2 
where £3(3;, ti, £2) and £4(2;, ii, £2) are any pair of functions identically satisfying the condition 

= -x + - \t\H + \t\ + hU 

solve both H51\) and {5. 



6 Conclusions 

In this article we presented a novel method of obtaining multicomponent Harry Dym hierarchy (both its 
local and nonlocal part) as well as wide classes of its solutions, from a family of finite dimensional separable 
systems (Stackel systems of Benenti type). This method has been previously applied to coupled Korteveg- 
de Vries hierarchy where it produced novel rational solutions and also a family of implicit solutions. In the 
case of cHD hierarchy discussed here, the method produces among others rational and implicit solutions 
in case of nonlocal hierarchy and implicit solutions of the local part. In addition, the method produces 
wide families of other solutions that are to be exploited elsewhere. It also indicates the existence of 
common solutions of local and nonlocal cHD systems. 

Our method can hopefully be extended to other systems, for example by taking more general separa- 
tion relations than relations 
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7 Appendix 



We prove here TheoremQjJ] We start with the case (3 — 0. For j3 = the solutions ([35]) are just solutions 
with our choice of m and fc, namely m — —N, k — 0. The functions 

qi(x = h,t2, . . . ,t n ), ■ ■ ■ ,q n (x = h,t2, ■ ■ ■ ,t n ) (58) 

obtained from (155]) (with (3 = 0) through ([7]) satisfy thus all n — 1 Killing systems (|T7| . Moreover they 
satisfy all the equations (|23l) and thus also all the equations ([25)1 used in our elimination procedure. This 
means that we are free to use any part of (|23l) or (1251) to perform an elimination of variables in (|17p . Such 
elimination thus leads to new equations that are satisfied by those functions from the set (|S"5|) that survive 
the elimination. Now, we know that replacing the variables gjv+i, ■ ■ ■ ,q n in the first N components of 
the first s — 1 = n — N equations (flT|) by the functions given by (f2"5j) leads to the first s — 1 flows of the 
hierarchy (|27|) . That means precisely that the first N functions in ([58)) 

qi{x = ti,t 2 , • • • ,t n ), . . .,q N (x = ti,ta, . . .,t n ) (59) 

satisfy the first s — 1 = n — N equations in (|27l) . We will now show that they actually solve the first n — 1 
equation in (|27p . Consider the next flow q ts+1 = Z^ +1 [q] in (l27l) . In order to obtain this flow, we have to 

perform the elimination of variables gjv+i, • ■ • , q-n, q-n+i in the flow qt 3+1 = Z™+1 [q] through (|25|) written 
for n + 1 instead of n i.e. obtained from solving (|24[) . This elimination is therefore performed with the 
help of the same functions qi — qi\q\. . . ■ , <zat] as for n plus a new function q n +i = f n ~N+2 [qi, ■ ■ ■ , Qn]' 
However, [q]) 3 = [q]) J for aU j = 1,...,N - 1 (it follows from $Tft) while (Z^l [q]) N 

contains the additional variable q n +i not present in (Z" +1 [q]) N ■ It means that solutions (|59[) will certainly 
satisfy the first N — 1 components in g ts+1 = Z^ +1 [q]. Further, since En+i (L" +1 ^ Ar '°) = —\q n +i.xx + 
Em (L n >~ N ^ , the function q n +\ — f n —N+2 fell • • • ■> Qn] is (after choosing both integration constants equal 
to zero) identically equal to zero on the solutions (|59p . That means that on the solutions (|5"9"]) we have 
(Z™+i [<?]) = {Z r s l + i [q]) N which means indeed that ([59^1 solves q ts+1 = Z^ +1 [q]. By expanding this 
argument, the functions q n +x , q n +2 ; • • • , q-n+N-i obtained from (J53J) with n replaced by n' = n + N — 1 
i.e. from the n' — N — n — 1 equations 

— \q n +N-l,xx + Vn-llqii —iQn+N-s] = 0) 

— \qn+N~2,xx + CPn-zlqii •••) 9«+JV-3] = 0, 

(60) 

~\qN+2,xx + ip2[qi,—,qN+i] = o 
\qN+i,xx + <Pi[qi, -,9jv] = 0. 

(which are necessary to obtain the first n—1 flows of ([27]) ) are identically zero on the solutions (|59|) which 
leads to the conclusion that ([59)) indeed solve the first n—1 equations of ([27]) . 

Assume finally that < /5 < n — 1. The functions (|35|) are then the complete solution (as usual, 
through the map ©) of all the Euler-Lagrange equations Ei{L n ^~ N ~^) associated with the Lagrangian 
£n,p-N,-p _ ^ g sucnj they still must solve all the Killing systems (|T7|) . However, since Ei(L n '@~~ N '~ l3 ) — 
Ei + p(L n ~ N '°) for i = 1, . . . , n — f) due to (|20]). for any /3 > 1 we lose the first ^ — 1 equations in (JHOJ) 
which means that our proof works only for the first n — (3 flows in (|27[) - we simply can not " blow up" n 
to n 1 = n + N - 1 but only to n" = n + N - 1 - (3. 
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E N+1 (L'^- N A = 
E N+2 (L™'.-^ ) = 

K>-i (L n '<- N <°) = 

E n , (L n '>- N A = - 
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